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Abstract 

The effect of local density of states (LDOS) fluctuations on the dynamics 

of a Kondo impurities in a small metallic point contact (PC) is studied. To 

estimate the spatial and energy dependent LDOS fluctuations we investigate 

a model PC by means of a transfer matrix formalism. For small PC's in the 

nanometer scale we find that near to the orifice strong LDOS fluctuations 

develop. These fluctuations may shift the Kondo temperature by several 

orders of magnitude, and result in a strong broadening of the PC Kondo peak 

in agreement with the results of recent measurements. 
PACS numbers:72.10.Fk, 83.30.H, 73.40.G 
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Magnetic impurities have been studied in the past few years very successfully using 
ballistic point contacts (PC's) Very recently a thorough reinvestigation of the 

Kondo effect has been carried out by Yanson et al. using small metallic PC's with the 
diameters in the mesoscopic range, d ~ 20A. Surprisingly, in contrast to measurements in 
thin metallic films and quantum wires, where a suppression of the Kondo effect has been 
found H, Yanson et al. observed a strong broadening of the zero bias Kondo peak which has 
been explained by the increase of the Kondo temperature of the impurities in the contact 
region p. 

In this Letter we report on the first study of the local density of states (LDOS) and 
their effect on the Kondo resonance in small ballistic PC's. We find surprisingly strong 
LDOS fluctuations in the contact region which are generated by scattering of the conduction 
electrons at the surface of the PC and give a natural explanation to the broadening of the 
Kondo peak in very small PC's ||. 

The effect of LDOS fluctuations on the Kondo resonance has already been studied in 
the case of disordered metals J7|. However, Dobrosavljevic et al. calculate the effect of the 
fluctuating LDOS only in the leading logarithmic order and their results are only valid in 
the case where the characteristic energy scale of the fluctuations, e c is of the same order of 
magnitude as the bandwidth cutoff, D ~ e c . As we shall see in a PC e c -C D thus we had to 
treat the fluctuations with more care. Furthermore, we also constructed the next to leading 
logarithmic scaling equations wich goes beyond the theory of Ref. 0. 

To study the effects of the LDOS fluctuations on a Kondo impurity we investigate the 
dynamics of a magnetic impurity in the contact region. For this purpose we use a simplified 
Hamiltonian. The first part of our Hamiltonian describes the conduction electrons: Hq = 
J2 n ,a e nO'ta a ncT, where the operator a+ CT (a na ) creates (annihilates) a conduction electron with 
spin a, wave function (p n , and energy e n . The states (p n are not momentum eigenstates 
but rather appropriately chosen scattering states of the PC. The interaction part of the 
Hamiltonian can be written as 
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flint = Yl <P*n( r ) J Vn>(r) a+ a ai a ,a n/(T/ S\ (1) 

where the impurity is at the position r, J is the strength of the local exchange interaction, 
and a and S denote the spin of the conduction electrons and the impurity, respectively. 

Far from the orifice the LDOS must approach its bulk value: £?(r, e) — > f?buik( e )- For 
the sake of simplicity we assume that the bulk DOS is constant between the high- and low- 
energy cutoffs, ^buik(e) = £?o- However, approaching the contact region g(r,e) deviates from 
its bulk value and random spatial and energy dependent LDOS fluctuations appear. 

To estimate the amplitude of these fluctuations we carried out numerical calculations 
for a free electron model PC consisting of two infinite half-spaces connected by a cylinder 
of radius R and length L |. To simplify the calculations we used angular momentum 
eigenstates propagating along the axis z of the PC of the form: 

V±,ex m (ri z, <p) = e^ m J m {\r)e ±lk ^ z , (2) 

where m is the angular momentum around the axis z and cylindrical coordinates have been 
used. The signs ± correspond to right- and left-going states, respectively, e denotes the 
energy of the conduction electron, and J m stands for the m'th Bessel function. The z 
component of the momentum in Eq.(^), k z , can be expressed by the energy e and the radial 
momentum A of the electron as k z = (2m e e — A 2 ) 1 / 2 and k z = i(X 2 — 2m e e) 1//2 for y / 2m e e\ and 
y / 2m e e < A, respectively ( m e denotes the electron mass). Since we require the vanishing of 
the wave functions at the boundaries, A is a continuous parameter in the infinite half-spaces 
while it takes discrete values inside the tube. 

The scattering matrices of the problem can be constructed by matching the wave func- 
tions and their derivatives at the two ends of the cylinder P,|T0fl. Having obtained the 
scattering matrices one can proceed by constructing all the scattering states and evaluating 
the different physical quantities like the LDOS or the conductance of the PC 0. 

Fig. |l].a shows the calculated LDOS fluctuations inside the tube for a PC with R = 15A 
and L = 15A at the point r = 7 A and z = 7 A, where the coordinate z is measured from 
the left wall of the PC. Each time a new conduction channel is opening a peak appears in 
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the function g(e, r, z). As one can see in Fig. |I|.b the strong interference effects cause strong 
fluctuations even for a fixed energy if the spatial coordinate is varied. 

In order to explore the effect of these LDOS fluctuations on the Kondo impurity we 
construct the next to leading logarithmic scaling equations for an arbitrary LDOS using the 
multiplicative renormalization group technique [pLl[J . Following Abrikosov |TJ|] we describe 



the impurity spin dynamics by means of a pseudofermion field, S l — > b+S l ss ,b s ', where the 
operators 6+ (b s ) create (annihilate) a pseudofermion corresponding to the spin state S z = s. 
The pseudofermion Green's function, Q(uj), and the vertex function, T nn /(uj) can be intro- 
duced in the usual way fl2fl . Due to the special structure of the interaction the self-energy 
of the conduction electrons [the vertex function] factorizes: E nn /(o;) = </?*(r)£(u;)</? n /(r) [ 
^nn'iyJi) = V 9 n( r )r(^j) ( / 9 r i '( r )]) where and T(ui) depend only on the LDOS at the posi- 
tion of the impurity. Then the multiplicative renormalization group equations can be written 



in the form [11 



g(u,f,D') = Z(j,D/D') G(u,j,D), 

r(u i7 f,D') = z^d/d')- 1 r(ui,j,D) , (3) 

where D' is the scaled bandwidth, Z is the pseudofermion wave function renormalization 
factor, and f denotes the dimensionless scaled coupling, f = J(D')gb ll i] 5 _(ep), ep being the 
Fermi energy. 

Similarly to the case of Kondo impurities or two level systems in a smooth LDOS 
|TT| , p!3|JT7| the leading logarithmic terms in the scaling equation arise from the second order 
vertex corrections while the next to leading logarithmic terms are generated by the second 
order pseudofermion self energy corrections and the third order vertex corrections. Having 
calculated these diagrams one can easily generate the scaling equations using Eq. (|) and 
after some lengthy algebra one obtains: 



(1 + £AD) 2 

(R(OR(-D)+R(-OR(D)) , (4) 



where R(£) = £>(£)/£buik(e.F), £ being the energy of an electron measured from the Fermi 
energy, £ = e — ep and we introduced the scaling variable x = ]n(D /D), D w being 
the initial bandwidth cutoff. Naturally, for a system without LDOS fluctuations R = 1 and 
Eq. (f|) reduces to the usual next to leading logarithmic scaling equation |TT| . 



To see how the LDOS fluctuations modify the Kondo temperature we drop the next to 
leading logarithmic term in Eq. (f|) and calculate the leading logarithmic Kondo temperature, 
which is associated to the divergence of the scaled coupling: j(D = Tk) = oo. Assuming 
that both the bulk and the fluctuation induced Kondo temperatures, T K and Tk are small 
enough to satisfy g(£ = Tk) ~ f?(£ = T K ) w g(ep) one easily obtains the following expression 
for the ratio of the two Kondo temperatures 

t\ = exp { £ % {6m + sr (-^} > ( s ) 

where 5R(£) = R(£) — RbuikiQ- The important feature of Formula (pD is the appearance 
of the weight l/£ in the exponent. The physical meaning of this factor is that the Kondo 
resonance is mainly formed by the low-energy electron-hole excitations. The contribution 
of high-energy electron-hole excitations is also significant but it is suppressed compared to 
that of the low-energy ones. Roughly speaking most of the sites with Sg(ep)0 will have an 
increased Tk while impurities for which Sg(ep) < will tend to have a Tk decreased. 

The relative Kondo temperature Tk/T k can be easily estimated by substituting the ratio 
5R(£) of the free electron calculations into Eq.(^[). The calculated average ratio < T K /T K 
as a function of the radius R of the PC for fixed length L = 5 A is shown in Fig. 2(a). The 
increased average Kondo temperature, < Tk ~ 1 — 10K, is orders of magnitude larger than 
bulk Kondo temperature, T K = O.Oli^, and is in reasonable agreement with the experiments 
(see the diamonds in the Figure) [Tj|. A rough fitting of our data with the experimentally 
found power law dependence < Tk ~ d~ a gives an exponent a = 2.2 ± 0.5 which agrees 
qualitatively with the experimental exponent, a = 2. (Unfortunately, the available PC 
diameters were limited by our computer capacity and only a rough estimation of a could 
have been carried out.) As a comparison we also show the average < Tk/T k calculated 



by integrating the next to leading logarithmic scaling equations, Eq.(f|). As one can see in 
Fig. 2(a) there is no significant difference between the next to leading logarithmic and the 
leading logarithmic results. 

Determining the resistivity contribution of a magnetic impurity at a site r in an ultrasmall 
PC is a very complicated task because at the length scales involved in the problem the 
usual Boltzman equation description [|14]] breaks down, and presently there is no theory 
available which could take into account both the geometrical effects (i.e., the strong quantum 
interference) and the strongly correlated behavior of a Kondo impurity in the PC. Therefore 
we proceed in a semiquantitative way and we estimate the scattering rate 1/t(u) of a 
conduction electron passing through the orifice. This is connected to the imaginary part 
of the conduction electrons' self energy, 1/t(u) = —2 Im £(&>), and can be measured 
directly in a PC experiment |||14]. Since the conduction electrons' Green's function remains 



unrenormalized in the dilute impurity limit this scattering rate is simply proportional to 
1/t(uj) ~ 2nS(S + 1)(q(€f)J(uj)) 2 6f, where J(u) denotes the scaled coupling at the energy 
D = uj. As in the unitary limit the contribution of an impurity in the contact region to 
the conductance of the PC is approximately ~ —e 2 /h in the u = eV — > limit , and 
q(ef)J(u) — > 1 for uj — > (see Eq.(fJJ)), the change in the conductance of the PC due to 
magnetic scattering can be estimated as 



e 2 

AG(eV) = -—n cg 2 (e F ) < J 2 (max{eV, T}) , (6) 

where c is the concentration of the impurities, < ... denotes the average over the contact 
region and Q = 8R 3 /3 is the effective volume of the PC ||. 

The calculated amplitudes of the impurity contributions to the zero voltage conductance 
of the PC as a function of the system size are compared to the experimental data in Fig. 2(b). 
The average in Eq. ([]) has been carried out over 40 randomly chosen impurity positions. 
Both the amplitude of AG(R) and its sample to sample fluctuations are in very good agree- 
ment with the experiments. We stress at this point that there is no free parameter in Eq. (^|) 
exept for the length of the PC which hardly influenced our results. To show that there is a 
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striking size effect in Fig. 2(b) we also plotted the calculated amplitude of the Kondo signal 
with impurities having the bulk Kondo temperature. Clearly, both the amplitude and the 
size dependence (AG ~ d 2 ' 2 ) of the experimental and our calculated fluctuation dependent 
Kondo conductances are quite different from the one we obtained by assuming the bulk 
Kondo temperature (AG ~ d 3 ). 

The size dependence of the Kondo resonance can be understood as follows. Decreasing 
the contact size strong LDOS fluctuations evolve in the contact region. The smaller the 
contact size the larger the fluctuations become, and thus for small contact sizes the measured 
few impurities (i.e., the impurities in the contact region) have a very broad distribution of 
Kondo temperatures. Since in the temperature range ~ IK where the measurements were 
performed || magnetic impurities with a small Kondo temperature, T% pa T£ ~ O.Oli^, 
give no significant contribution to Eq.(^jj) the differential conductance is always dominated 
by the impurities with the largest Kondo temperatures. These arguments suggest that the 
sum in Eq.(|6|) is dominated by the magnetic impurities having large Kondo temperatures. 
Therefore our theory gives a natural explanation why the applied magnetic field could not 
destroy the zero bias Kondo anomaly for small samples ||. 

One can also show analyzing formula (^) that for a fixed PC shape the relative increase 
of the Kondo temperature, Tjc/T^, depends essentially on the bulk Kondo temperature T£, 
and for alloys with very small T£ the ratio T^/T£ can be much larger than for alloys with 
large T%. This explains why no significant broadening of the Kondo peak has been observed 
in CuFe PC's (T£ pa 20K) [|J, while an enormous broadening was found in AuMn and 
CuMn {T* K < OMR) H§. 

According to the picture proposed in this Letter one expects that any physical quantity 
which is sensitive on the LDOS will have a distribution in the contact region due to the 
presence of strong fluctuations. As an example we mention the PC spectrum of fast two 
level systems (TLS) [|H|[nj. Since the dimensionless TLS - conduction electron couplings 
are sensitive on the LDOS [[T^,|T^] and Tg is in the range of ~ IK |I5| we expect a slight 



broadening of the zero bias TLS peak, which has indeed been observed in metallic break 



junctions |L6] . 

In our simple model the LDOS fluctuations were basically connected to the openings 
of new conductance channels through the PC [[J. Such LDOS fluctuations, however, may 
also be generated by random scattering at the boundary of the PC. Effectively, it has 
been found recently that random scattering at the surface of a PC is able to produce huge 



LDOS fluctuations in the contact region |20[. Therefore we think that our conclusions are 
qualitatively independent of the special model considered. 

The case of ultrasmall PC's with magnetic impurities f| should be contrasted to the 
experiments performed on thin metallic films and wires |J. In the latter case only alloys with 
high Kondo temperatures (CuFe, AuFe) were investigated, where the LDOS fluctuations 
do not modify the Kondo temperature very much. Instead, it seems to be that spin orbit 
interaction induced spin anisotropy develops which leads to the blocking of the spin flip 



processes |2T| and a suppression of the Kondo signal. On the contrary for thin AuMn films, 
where T£ is very small, we expect that the Kondo temperature can be increased due to the 
surface induced LDOS fluctuations by several orders of magnitude similarly to the case of 
PC's. 

In summary, we have proposed that the recently measured anomalous Kondo tempera- 
tures in small CuMn PC's are due to the strong LDOS fluctuations in the contact region 
generated by the scattering of the conduction electrons at the surface of the PC. Carrying 
out a model calculation we found that these fluctuations are really able to increase the 
Kondo temperature by several orders of magnitude in agreement with the measurements. 
Our results also suggest that the Boltzman equation approach becomes inadequate for ul- 
trasmall PC spectrometry, where quantum interference and fluctuation effects may become 
extremely important. 

The authors are deeply indebted to A. Zawadowski, I. K. Yanson, Jan van Ruitenbeek, 
N van der Post, and I. V. Kozub for helpful discussions. This research has been supported 
by the Hungarian Grants OTKA 7283/93 and OTKA F016604 . 
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FIGURES 

FIG. 1. (a): Fluctuations of the LDOS g(e,z,r) inside the tube for a PC with R = lhA and 
L = IbA at the point r = 7 A and z = 7 A. The fluctuating LDOS does not integrate exactly to 
its bulk value (dashed line) because the hard wall of the PC pushes the conduction electrons in 
the inside region of the tube, (b): g(e,z,r) for the same PC at z = 7 A and energy e = 7eV as a 
function of the radius r. 

FIG. 2. (a) The average relative Kondo temperatures < T^/T^ as a function of the radius R of 
a CuMn PC with L = 5A. The cutoff and the bare couplings have been chosen to be Dq = Q.SeV 
and jo = 0.032. Crosses and diamonds denote the results obtained in the leading and in the next 
to leading logarithmic orders, while data points indicated by boxes have been obtained from the 
results of Ref . || . The large fluctuations are due to the sensitivity of the interference pattern to the 
geometry of the PC. (b) Size dependence of the amplitude of the dimensionless Kondo conductance 
Ag = AG h/e 2 for the same contact. Diamonds denote the experimental data taken from Ref. || 
while our results are indicated by crosses. The impurity concentration and the temperature were 
c = 0.1% and T = 0.05K, respectively. The dashed line indicates the results without LDOS 
fluctuations (g ~ R 3 ) while the continuous line corresponds to the best fit to the data of Ref. ||: 
g ~ i? 217 . 
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FIGURES 
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FIG. 1. G. Zarand and L. Udvardi: "Enhancement of the Kondo temperature of magnetic 
impurities in metallic point contacts..." 
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FIG. 2. G. Zarand and L. Udvardi: "Enhancement of the Kondo temperature of magnetic 
impurities in metallic point contacts..." 
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